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Abstract 

The inclusive production of D* 1 * 1 mesons in two-photon collisions is measured with 
the ALEPH detector at e + e~ centre-of-mass energies from 183 GeV to 209 GeV. A 
total of 360 ± 27 D*^ meson events were observed from an integrated luminosity 
of 699pb~L Contributions from direct and single-resolved processes are separated 
using the ratio of the transverse momentum p]p* of the D* 1 * 1 to the visible invariant 
mass VLvis of the event. Differential cross sections of D*^ production as functions of 
pY* and the pseudorapidity \rp* | are measured in the range 2GeV/c < p{p* < 
12GeV/c and \rp \ < 1.5. They are compared to next-to-leading order (NLO) 
perturbative QCD calculations. The extrapolation of the integrated visible D*^ 
cross section to the total charm cross section, based on the Pythia Monte Carlo 
program, yields cr(e + e" -> e+e"cc) <v ^ >=197GeV = 731 ± 74 stat ± 47 syst ± 157 extr pb. 
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1 Introduction 

Heavy flavour production in two-photon events at LEP 2 centre-of-mass energies is dom- 
inated by charm production processes in which both of the photons couple directly (di- 
rect processes) or in which one photon couples directly and the other appears resolved 
(single-resolved processes) (Fig. [TJ [lj. These two contributions are of the same order 
of magnitude within the experimental acceptance. Because the single-resolved process 
is dominated by jg fusion, the measurement of the cross section can give access to the 
gluon content of the photon. Moreover, the large masses of the c and b quarks provide 
a cutoff for perturbative QCD calculations, allowing a good test of QCD predictions for 
the corresponding reactions. Contributions from processes in which both photons appear 
resolved (double-resolved processes) are suppressed by more than two orders of magnitude 
compared to the total cross section [T]. The production of b quark is expected to be 
suppressed by a large factor compared to charm quark because of the heavier mass and 
smaller absolute charge. 

In the present analysis charm production is measured in two steps. A high-purity 77 
sample is first selected, then examined for its charm content via reconstruction of D* + 
mesons in their decay to D°7r + . This letter is organized as follows. A short descrip- 
tion of the ALEPH detector is given in Section |2j Monte Carlo simulations for signal 
and background processes are described in Section El In Section HI event selection and 
reconstruction of D* + mesons are discussed. The results of the analysis are presented 
in Section |S| Finally, in Section |H1 a summary is given. Throughout this letter charge- 
conjugated particles and their decays are implicitly included. 

2 ALEPH Detector 

The ALEPH detector has been described in detail in (21 E] • Here, only the parts essential 
to the present analysis are covered briefly. The central part of the ALEPH detector is 
dedicated to the reconstruction of the trajectories of charged particles. The trajectory 
of a charged particle emerging from the interaction point is measured by a two-layer 
silicon strip vertex detector (VDET), a cylindrical drift chamber (ITC) and a large time 
projection chamber (TPC). The three tracking detectors are immersed in a 1.5 T axial 
magnetic field provided by a superconducting solenoidal coil. Together they measure 
charged particle transverse momenta with a resolution of Spt/pt = 6 x 10~ 4 p t ©0.005 (p t in 
GeV/c). The TPC also provides a measurement of the specific ionization dE / dx meas . An 
estimator Xh = (dE/dx meas — d£ , /dx e xp,h)/o"cxp,h is formed to test a particle hypothesis, 
where dE / dx exp ^ an d o~ exp ^ denote the expected specific ionization and the estimated 
uncertainty for the particle hypothesis h, respectively. A mass hypothesis may be tested 
by means of the Xh values themselves or by calculating Xh confidence levels P^. 

Photons are identified in the electromagnetic calorimeter (ECAL), situated between 
the TPC and the coil. The ECAL is a lead/proportional-tube sampling calorimeter 
segmented in 0.9° x 0.9° projective towers and read out in three sections in depth. It 
has a total thickness of 22 radiation lengths and yields a relative energy resolution of 
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0.18/v^E + 0.009, with E in GeV, for isolated photons. Electrons are identified by their 
transverse and longitudinal shower profiles in ECAL and their specific ionization in the 
TPC. 

The iron return yoke is instrumented with 23 layers of streamer tubes and forms the 
hadron calorimeter (HCAL). The latter provides a relative energy resolution of charged 
and neutral hadrons of 0.85 /y/E, with E in GeV. Muons are distinguished from hadrons 
by their characteristic pattern in HCAL and by the muon chambers, composed of two 
double-layers of streamer tubes outside HCAL. 

Two small-angle calorimeters, the luminosity calorimeter (LCAL) and the silicon lu- 
minosity calorimeter (SICAL), are particularly important for this analysis to veto events 
with detected scattered electrons. The LCAL is a lead/proportional-tube calorimeter, 
similar to ECAL, placed around the beam pipe at each end of the detector. It monitors 
angles from 45 to 160 mrad with an energy resolution of 0.15W^B(GeV). The SICAL 
uses 12 silicon/tungsten layers to sample showers. It is mounted around the beam pipe 
in front of the LCAL, covering angles from 34 to 58 mrad, with an energy resolution of 
0.225 v /E(GeV). 

The information from the tracking detectors and the calorimeters are combined in an 
energy-flow algorithm [3]. For each event, the algorithm provides a set of charged and 
neutral reconstructed particles, called energy-flow objects. 

3 Monte Carlo Simulations 

In order to simulate the process e + e~ — > e + e~77 — > e + e~cc — > e + e _ D* ± X, the leading- 
order (LO) PYTHIA 6.121 Monte Carlo [I] is used. Events are generated at e + e~ centre- 
of-mass energies ranging from 183 GeV to 209 GeV using the corresponding integrated 
luminosities for weighting. Two different samples, direct and single-resolved processes, 
were generated for each of the considered D* + decay modes using matrix elements for 
the massive charm quark. The charm quark mass m c is chosen to be 1.5 GeV/c 2 and the 
parameter Aqcd is set to 0.291 GeV/c 2 . The 77 invariant mass is required to be at 
least 3.875 GeV/c 2 , which is the DD threshold. In order to ensure that both photons are 
quasi-real, the maximum squared four-momentum transfer Q^ax * s limited to 4.5 GeV/c 2 . 
In the single-resolved process, the SaS-lD j3] parametrization is used for the partonic 
distribution of the resolved photon. The Peterson et al. parametrization [6J is adopted as 
the fragmentation function of the charm quark with the nonperturbative parameter e c = 
0.031. The background process e + e~ — > e + e~77 — > e + e~bb is simulated using PYTHIA 
6.121 with Wj-y being required to be at least 10.5GeV/c 2 , which is the BB threshold. 
The b quark mass is set to 4.5GeV/c 2 . Again the Peterson et al. parametrization is 
adopted with eb = 0.0035. Other possible background processes have been simulated 
using appropriate Monte Carlo generators as listed in Table [T] 
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4 Event Selection and Reconstruction of D* + Mesons 
4.1 Selection of 77 Events 

The data analyzed were collected by the ALEPH detector at e + e~ centre-of-mass energies 
ranging from 183 GeV to 209 GeV with an integrated luminosity C = 699 ptr 1 . The event 
variables used for the event preselection are based on the ALEPH energy-flow objects. 
The following cuts, derived from Monte Carlo studies, were applied to select two-photon 
events. 

• The event must contain at least 3 charged particles. This cut reduces the background 
from leptonic events. 

• The visible invariant mass W w - ls of the event must lie between 4 GeV/c 2 and 55 GeV/c 2 
while the total energy of charged particles E c ^ should not exceed 35 GeV in order 
to reject the e + e~ annihilation background. 

• The visible transverse momentum p t vis of the event is required to be less than 
8GeV/c, as the p t)Vis distribution has a much longer tail for all considered back- 
ground processes. 

• To reject further background processes a cut combining the number of charged tracks 
and the visible energy E v i S of the event is applied: < 40 — |£' v i s (GeV). 

• Finally, in order to retain only events with almost on-shell photons an anti-tagging 
condition was applied, i.e., tagged events were rejected. A tag in this analysis is 
defined as an energy-flow object in the luminosity calorimeters (LCAL and SICAL) 
with an energy of at least 30 GeV. 

This selection retains a sample of 4.9 million events. Monte Carlo studies of possible 
background sources predict a 77 purity of 98.8%. 



4.2 Reconstruction of D* + Mesons 

Charm quarks are detected using exclusively reconstructed D* + mesons which decay via 
D* + — > D°7r + , with the D° being identified in three decay modes, (1) K~7r + , (2) K~7r + 7r°, 
and (3) K~7r + 7r~7r + . As a basis for possible and 7r =l= candidates reconstructed tracks 
of charged particles which fulfill the following quality conditions are used: 
p >100MeV/c (momentum of track), 

I do I < 2 cm (distance to beam axis at closest approach), 

\z 1 <8cm (z coordinate at closest approach), 

N TPC > 4 (number of hits in TPC), 

|cos#| < 0.94 {0 = polar angle with respect to beam axis). 

A track surviving these cuts is classified as a kaon if the measured specific energy loss 
dE/dx of the track is consistent with the expectation value for the kaon mass hypothesis, 
i.e., if the corresponding confidence level Pk is greater than 10%. The track is classified 
as a pion if P n is at least 1%. Thus, each track can be flagged as a kaon or pion or both 
or neither. 
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The 7T° candidates are formed from pairs of photons found in ECAL with an energy of 
at least 250 MeV each and an invariant mass within 85MeV/c 2 of the nominal 7r° mass. 
In order to improve the energy resolution of these 7r°'s the energies of the photons are 
refitted using the 7r° mass as constraint. If the confidence level of this fit is greater than 
5% and if | cos^ol < 0.93, where 9 n o is the polar angle of the n° candidate with respect 
to the beam axis, the 7r° candidate is retained. 

The D° candidates are formed from appropriate combinations of identified kaons and 
pions according to three considered decay modes. The D° candidate is retained if it has 
an invariant mass within 20MeV/c 2 , 65MeV/c 2 , and 20MeV/c 2 of the nominal D° mass 
for decay mode (1), (2), and (3), respectively. These mass ranges correspond to about 
three times the mass resolution. In order to reduce the combinatorial background in mode 
(3), the four tracks composing the D° are fitted to a common vertex and the confidence 
level of this fit is required to be greater than 0.2%. The combination of each D° with one 
of the remaining ir + candidates is considered to be a D* + candidate. In order to reduce 
combinatorial background from soft processes and to limit the kinematic range of the D* + 
to the acceptance range of the detector with reasonable efficiency, cuts were applied to 
the transverse momentum p t and the pseudorapidity rj = — ln(tan(6 l /2)) of the D* + : 

2GeV/c<p t D * + < 12GeV/c, |r/ D * + | < 1.5 . (1) 

If there are several D* + candidates found in one event the corresponding D° candidates 
are compared in mass and only the candidate with D° mass nearest the nominal D° mass 
is retained. If two or more D* + candidates share the same D° candidate, all of them are 
retained. Figure 121 shows the mass difference Am = m^*+ — m-^o for the selected D* + 
candidates for all three decay modes together. The spectrum rises at the lower threshold 
given by the pion mass. A clear peak is seen around 145.5 MeV/c 2 . In order to extract the 
number of D* + events the data distribution is fitted with the following parametrization: 



F(Am) = N 



1 f I f Am- 145.5 MeV/c 2X 2 



2ira 



exp< > + C (Am — m 



• (2) 



In order to exclude systematic binning effects an unbinned maximum likelihood fit is 
performed where C and P are used as free parameters. The normalization N follows from 
the constraint that the integral of F(Am) over the range of the fit, 130MeV/c 2 < Am < 
200 MeV/c 2 , must be equal to the number of entries in the histogram. The width a of the 
Gaussian describing the peak is fixed to 0.5 MeV/c 2 , as determined in Monte Carlo. The 
number of D* + events is then obtained by integrating the Gaussian part of @2J) in the range 
of 145.5 MeV/c 2 ± 3a. As the result a total of 360.0 ± 27.0 stat D*+ events are observed 
for all three D* + decay modes together. Among the possible background processes, only 
the contribution from 77 — > bb — > D* ± X is found to be sizeable. This contribution is 
estimated to be 20.5 ± 1.6 s tat D* + events from a 77 — > bb — > D* ± X Monte Carlo sample 
and the total cross section cx(e + e~ — > e + e~bb) measured in [TT]. After subtraction of this 
background, a total of 339.5 ± 27.0 stat D* + events are found in the data sample analyzed. 
The mass difference distributions for three channels separately are shown in Fig. OJ 
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5 Cross Section Measurements 



5.1 Relative Fractions of Direct and Single-resolved Contribu- 
tions 

As mentioned in the introduction, open charm production in 77 collisions is dominated 
by contributions from direct and single-resolved processes. In the direct case the cc pair 
makes up the final state of the 77 system (in LO) whereas in the single-resolved case the 
partons of the resolved photon (photon residue) in addition to the cc pair make up the final 
state. The transverse momentum p®* + of the D* + is correlated with the invariant mass 
of the cc system and the total visible invariant mass W v i S is in turn correlated with the 
invariant mass of the total 77 system. The ratio p®* + /W v[s should therefore be distributed 
at higher values for the direct case compared to the distribution of single-resolved events. 

FigurelUshows the distribution of p^* + /W Y - ls in data for all events found in the signal re- 
gion of the mass- difference spectrum. Combinatorial background has been subtracted us- 
ing events of the upper sideband 0.16 GeV/c 2 < Am < 0.2 GeV/c 2 of the mass-difference 
spectrum. Background from bb production has also been subtracted. The relative frac- 
tions are determined by fitting the sum of the direct and single-resolved Monte Carlo 
distributions to data with the relative fraction as a free parameter of the fit. The total 
number of entries in this Monte Carlo sum is required to be equal to the number of entries 
in the data distribution. The fit yields a direct contribution of r^r = (62.6 ± 4.2)% and 
a single-resolved contribution of r res = 1 — r dir = (37.4 ± 4.2)%. 



5.2 Differential Cross Sections 

Two differential cross sections for the production of D* + mesons are determined: the first 
one as a function of the transverse D* + momentum and the second as a function 

of pseudorapidity |t7 D + |. Both are restricted to the range defined in Eq. (0). The former 
is measured in three pf* + bins: [2-3], [3-5], [5-12] GeV/c, and the latter in three |?7 D * + | 
bins: [0-0.5], [0.5-1.0], [1.0-1.5]. All considered D* + decay modes were treated separately. 

The average differential cross section da/dp®* + for a given p^* + bin and |^ D * + | < 1.5 
is obtained by 



dp t D * + Ap»* + £B*B e pD * + 

ft 



(3) 



Analogously one obtains dcx/d|?7 D * + | for a given bin in |?7 D * + | and 2GeV/c < p^* + < 
12GeV/c 



where 

• -^found i s the number of D* + found in the considered bin after subtracting the bb 
background (determined as described in Section 14 .2|) with the width of the fitted 
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Gaussian (J2J) being fixed to 0.5 MeV/c 2 for decay modes (1) and (3) and 0.7MeV/c 2 
for decay mode (2), 

• Ap® , A|t] d * + | are the considered intervals in p®* + and |// D * + |, 

• £ = 699 pb -1 is the integrated luminosity of the data analyzed, 

• is the branching ratio BR(D*+ -> D tt + ) = (68.3 ± 1.4)% PH, 

• Bq is the branching ratio of the considered D° decay mode [THj . 

• e D »+ (e, D .+ ,) is the efficiency of reconstructing a D* + candidate in the given pP* + 

(|?7 D * |) bin in the considered decay mode. Since efficiencies are determined sep- 
arately for direct and single-resolved processes (e d n*+ and e r %+, respectively) the 

Pt Pt 

total efficiency is a weighted combination using the fractions as determined in Sec- 
tion EU 

£ D*+ = r dir£ d*+ + r res£ d*+ j 
Pt P t Pt 

, _ „ ,dir I res 

fc |jjD*+| — 'dirt|^D*+| i 'restjjjD'+i 

Tables El and H2 show the number of D* + mesons found in the chosen p^* + and |?7 D * + | bins 
respectively, as well as the derived differential cross sections da/dp^* + and do"/d|r] D * + | 
with their statistical and systematic errors. The resulting cross sections for the different 
D* + decay modes are consistent with each other for all bins in pf* + as well as in |// D * + |, 
taking into account the statistical uncertainties. The weighted average over all of the 
considered D* + decay modes is given in Table |U for each p^* + and |^ D * + | bin, where only 
the dominating statistical uncertainties are used for weighting. 



5.2.1 Systematic Errors on Differential Cross Sections 

The study of systematic errors was performed separately for each pf* + and |^ D * + | bin and 
for each of the considered D* + decay modes, unless otherwise specified. 

The systematic error introduced by the event selection was studied by varying the cuts 
within the resolution obtained from the Monte Carlo detector simulation. The systematic 
uncertainty was estimated from the resulting relative variation of the efficiency. This 
yields an uncertainty of 0.6%-6.4%, depending on the considered pf* + or |?7 D * + | bin and 
on the D* + decay mode. 

The selection of pion and kaon candidates depends essentially on the dE/dx measure- 
ment as well as on the expectation values dE j da;exp,h used to calculate the probability 
for a given mass hypothesis m^. The uncertainty of the dE/dx calibration changes the 
efficiency by 0.5%-5.8%. These deviations are used as an estimate of the systematic error. 

The systematic error due to the accepted mass range used to classify D° candidates 
was examined by comparison of the mass distributions of D° candidates which contributed 
to the D* + signal in data and Monte Carlo for each D° decay mode separately. A Gaussian 
fit was applied to these distributions. The fraction of the fitted Gaussian which lies within 
the accepted mass range differs between data and Monte Carlo by less than 0.6%. Thus, 
the uncertainty due to this source can be neglected. 

In order to estimate the error introduced by the method for extracting the number 
of D* + events (Section 14 .2j) the mean of the fitted Gaussian in Eq. (J2J) was varied by 
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±0.05MeV/c 2 , and the width was varied by 10% about its values as obtained in Monte 
Carlo. The resulting relative error on the efficiencies was 0.8%-2.1%. 

A variation of the interval that defines the upper sideband yields a variation in r^ iT 
of less than 0.05%. Hence, this source is negligible. The present analysis assumes the 
fraction rdi r to be constant over the considered kinematic range. Monte Carlo studies 
show a variation of this fraction of up to 12% in this range, depending on the bin in 
p®* + and |?7 D * + |. A relative uncertainty of 10% is therefore added in quadrature to the 
statistical uncertainty of r ros . A variation of r dir / ros within these uncertainties yields a 
variation in the cross section of 0.3%-3.4%, which is used to estimate the introduced 
uncertainty. 

The statistical error of bb background subtraction and the uncertainties of the total 
cross section cr(e + e~ — > e + e~bb) yield a systematic error of 1.2%-3.4% on the differential 
cross sections. 

The overall trigger efficiency of the selected D* + events is estimated to be consistent 
with 100% with a statistical uncertainty of 1%. Thus no correction is made for this source. 

The relative errors on the branching ratios given in ^0] are used to estimate the 
corresponding relative systematic uncertainties in the cross sections. 

Similarly the relative uncertainties in the efficiencies due to finite statistics in the 
Monte Carlo samples, 0.5%-2.3%, are taken into account. 

All systematic errors are assumed to be uncorrelated and therefore added in quadra- 
ture. Table El shows a summary of the systematic uncertainties. 

5.2.2 Comparison to Theory 

Figures El and El show the measured da/dp^* + and dcr/d|?7 D * + | in comparison to two dif- 
ferent NLO perturbative QCD calculations, the fixed-order (FO) NLO (also known as 
massive approach) [T2] and the resummed (RES) NLO (massless approach) [T3|. In both 
cases, the charm quark mass m c is set to 1.5GeV/c 2 , the renormalization scale /iR and 
the factorization scale /ip are chosen such that /i| = = m\ = m 2 + p t (c) 2 , where 
p t (c) is the transverse momentum of the charm quark. For the resolved contribution the 
photonic parton densities of the GRS-HO parametrization are chosen ^3] in the FO NLO 
calculation, whereas the RES NLO uses GRV-HO [T^j. The fragmentation of the charm 
quark to the D* + is modelled by the fragmentation function suggested by Peterson et 
al. with e c = 0.035 in the case of FO NLO. The RES NLO calculation uses e c = 0.185, 
which was determined by using nonperturbative fragmentation functions fitted ^3] to 
ALEPH measurements of inclusive D* + production in e + e~ annihilation ^B] . The results 
of the two NLO QCD calculations are represented by the dashed lines (for RES NLO) 
and solid lines (for FO NLO) in both Fig. El and Fig. El In order to estimate the the- 
oretical uncertainties, the FO NLO calculation was repeated with the charm mass and 
the renormalization scale varied as described in the figures. The RES NLO calculation is 
also repeated using the AFG ^7j ansatz as an alternative for parton density function and 
varying the renormalization and factorization scales. The resulting theoretical uncertain- 
ties are indicated by the bands around the corresponding default values in Figs. El and 

El 

Altogether, the measurement of da/dpf* + seems to favour a harder p®* + spectrum 
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than predicted. The RES NLO calculation clearly overestimates the measurement in 
the low pf* + region, while the FO NLO calculation slightly underestimates it in the 
pf + > 3.0GeV/c region. The measured da/d|?7 D * + | is consistent with the almost flat 
distribution predicted by both NLO calculations, but the measurement of do"/d|r7 D * + | is 
again overestimated by the RES NLO calculation and somewhat underestimated by the 
FO NLO calculation. 

5.3 Visible Cross Section 

The visible cross section cr^* + (e + e~ — > e + e~D* + X) is calculated separately in the accep- 
tance range [Eq. ([JJ] for the three considered decay modes by 

°$? (e + e- - e + e"D*+X) = , (4) 

where the notation is as the same as in Eq. (J3J). The numbers of D* + found and the 
efficiencies of reconstructing a D* + candidate for direct and single-resolved processes are 
listed in Table |B1 together with the derived visible cross sections cr°* + (e + e~ — * e + e~D* + X) 
and their uncertainties for the three decay modes. The systematic error is determined 
in the same way as for differential cross sections ( Section 15.2.1)1 . The weighted average 
over all of the considered decay modes using the dominating statistical uncertainties for 
weighting is 

*3T(e + e- -> e+e-D*+X) = 23.39 ± 1.64 stat ± 1.52 syst pb . (5) 
The theoretically predicted cross section [12] is 

o3T (e + e- - e + e"D* + X) = 17.3+^ pb , (6) 
and is consistent with this measurement within the given uncertainties. 

5.4 Total Cross Section 

The total cross section for the reaction e + e~ — ► e + e~cc is given by 

a(e + e ^e + e cc) = -^ (r d irRdh + r res R res ) , (7) 

where the symbols are as follows: 

• o"°s + is the visible inclusive D* + cross section determined in the previous section; 

• P c -,d*+ is the probability for a charm quark to fragment into a D* + meson (taking 
the combined quantity P C ^ D *+ x BR(D*+ -> D tt + ) = 0.1631 ± 0.0050 from 18 and 
using BR(D*+ -> D tt + ) = (68.3 ± 1.4)% [ID] yields P C ^ D *+ = 0.2388 ± 0.0088); 

• the factor 2 in the denominator takes into account that, for the single inclusive cross 
sections, both the D* + and the D*~ mesons were counted; 
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• rd ir and r res are the fractions of the direct and single-resolved contributions in the 
considered acceptance range, as described in Section I3TT1 

• Rdir is the ratio 

D*+ 

T-) _ tot,dir 
u vis,dir 

of the total D* + cross section to the visible cross section in the range of Eq. for 
direct processes. It describes the extrapolation of the measured cross section to the 
total phase space available. R res is the corresponding quantity for the single-resolved 
case. 

Separate Monte Carlo samples are used to estimate -Rdir and -R re s for direct and single- 
resolved processes. The parameters used to determine -Rdir and -R re s are described in 
Section El This yields R dir = 12.74 ± 0.45 stat and R Tes = 18.62 ± 0.80 stat . 

The main theoretical uncertainties entering the calculation of the extrapolation factors 
stem from the uncertainty of the charm quark mass. A variation of the charm mass to 
m c = 1.3 GeV and m c = 1.7 GeV yields relative errors on R dir of ±10% and on i? res of 
+43% and —19%, respectively. 

In the single-resolved case an additional uncertainty enters -R res by the choice of the 
parton density functions describing the resolved photon. Alternatively to the default 
choice the GRV-LO parametrization [T5] was used to calculate R Tes . This yields a relative 
deviation of 12% and is added in quadrature to the other systematic uncertainties on R Tes . 
The following values are therefore obtained: 

.Rdir = 12.7 ± 1.3 

-Rres = 18-6^4 g ■ 

The uncertainties in rdir, c"vis + > and P c ^d*+, which are assumed to be uncorrelated, are 
taken into account in the estimation of the statistical and systematic error on the total 
cross section by Gaussian error propagation. This procedure yields a total cross section 
for the reaction e + e~ — > e + e~cc at e + e~ centre-of-mass energies yfs = (183 — 209) GeV, 
corresponding to the luminosity weighted average of 197 GeV, 

+ 157 

cr(e + e" -> e + e"cc) <v /5 >=197GeV = 731 ± 74 stat ± 47 syst _ 86cxtr pb . (8) 

Alternatively, the total cross section is determined by means of the NLO calculation 
referenced in the previous section; in this case the cross section is given by 

g (e+e-->e + e-cc)= ffi R M . (9) 

■^-^ c— >D*+ 



The value R to t = 22.2 is extracted from [12] by determining the ratio of the calculated 
total charm cross section to the charm cross section calculated for the visible D* + range 
considered in the present analysis. Variation of the parameters entering the calculation 
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yields deviations in the range from —33% to +72%, which are used as an estimate of the 
systematic error due to the extrapolation. This results in a total cross section 

cr(e + e" -> e + e"cc) <v ^ >= i 97GcV = 1087 ± 86 stat ± 70 syst ^ 3 gy pb . (10) 

The measured total cross section [Eq. (jSj)] is shown in Fig. [7| in comparison to the 
NLO QCD prediction of Drees et al. P and to results from other experiments |2"Ull21[l2*2"l 
122] • Within the uncertainties, this NLO QCD prediction is in good agreement with our 
measurement and others |2~4] . 

6 Conclusions 

The inclusive production of D* + mesons in two-photon collisions was measured using the 
ALEPH detector at LEP 2 energies in the reaction D*+ -> D°vr + . The D° mesons were 
identified in the decay modes K~ir + , K _ 7r + 7r°, and K~7r + 7r~7r + . A total of 339.5 ± 27.0 
D* + events from 77 — > cc was found in the kinematic region 2 GeV/c < p^* + < 12 GeV/c 
and |?7 D * + | < 1.5. 

The fractions of the main contributing processes, direct and single-resolved, were de- 
termined using the event variable pf* /W Y - IS to be = (62.6±4.2)% and r res = 1 — rdir = 
(37.4 ± 4.2)%, within the acceptance. 

The differential cross sections da/dp^* + and da/d|?7 D * + | were measured and compared 
to the fixed-order (FO) NLO QCD calculation |T2j and the resummed (RES) NLO QCD 
calculation ^3]. While the data show a slightly harder spectrum in the p^* + distribu- 
tion compared to both calculations, the almost flat distribution of dcr/d|?7 D * + | which is 
predicted by the NLO calculations for the visible D* + region is in agreement with the 
measurement. Overall, the measurements of da/dpf* + and dcr/d|?7 D * + | were slightly un- 
derestimated by the FO NLO calculation and overestimated by the RES NLO calculation. 

For the integrated visible D* + cross section a value of <x°* + = 23.39 ± 1.64 sta t ± 
1.52 syst pb is obtained which is consistent with the FO NLO calculation. 

The extrapolation of the visible D* + cross section to the total cross section of charm 

production introduces large theoretical uncertainties and has a large relative uncertainty. 

Using the LO calculation of the Pythia Monte Carlo we obtain 

_ _ +157 
cx(e + e" -> e + e"cc)< / i>=197GcV = 731 ± 74 stat ± 47 syst pb . 

A different method using the results from the FO NLO calculation yields a higher 
cross section and a larger error. 
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Table 1: Considered background processes and associated Monte Carlo generators 



Process 


Monte Carlo Generator 


e + e _ — > 


qq 


PYTHIA 5.7 m 


e + e~ — > 


T + T~ 


KORALZ 4.2 


e + e~ — > 


e + e~r + T~ 


PHOT02 |E] 


e + e~ — > 


w+w- 


KORALW 1.21 



Table 2: The numbers of D* + mesons found with \r] \ < 1.5 in bins of pjp for the 
three decay modes after background subtraction. The efficiency is listed separately for 
direct and single-resolved processes. The differential cross section in bins of the transverse 
momentum p^* + of the D* + for each considered D* + decay mode is given together with 
statistical and systematic errors. 



D*+ 

range 


iv found 


[GeV/c] 


D*+ -> (K-7r+)vr+ 


D*+ -» (K-7r+7r°)7r + 


D*+ -» (K~7r + 7r^7r+)vr+ 


2-3 


69.8 ±10.7 


18.7 ±6.2 


54.5 ±10.3 


3-5 


72.2 ±8.1 


29.0 ±7.8 


44.9 ±9.7 


5-12 


15.1 ±3.0 


20.9 ±5.7 


29.2 ±6.8 





Efficiency for direct process e d Q»+(%) 

pi. 


D*+ -> (K-7r + )vr+ 


D*+ -» (K-^+7r°)7r + 


D*+ -» (K~7r+7r^7r+)vr+ 


2- 3 

3- 5 
5-12 


27.96 ±0.13 
46.94 ±0.20 
48.73 ±0.34 


2.27 ±0.04 
6.83 ±0.10 
12.32 ±0.23 


11.66 ±0.09 
24.16 ±0.17 
30.13 ±0.33 








Efficiency for single-resolved process e T %+ (%) 

p± 


D*+ -> (K-7r + )vr+ 


D*+ -» (K-^+7r u )7r+ 


D*+ -» (K-7r + 7r-7r+)vr+ 


2- 3 

3- 5 
5-12 


26.81 ±0.12 
41.95 ±0.21 
34.59 ±0.41 


2.12 ±0.04 
6.17±0.10 
8.8 ±0.24 


10.49 ±0.09 
20.78 ±0.18 
19.83 ±0.36 








da/dpf + (pb/GeV/c) 


D*+ -> (K-7r + )7r+ 


D*+ -► (K-^+^>+ 


D*+ (K-7r+7r-7r + )vr + 


2- 3 

3- 5 
5-12 


13.80 ±2.12 ±1.04 
4.36 ±0.49 ±0.22 
0.27±0.05±0.01 


12.70 ±4.21 ±1.20 
3.32 ±0.90 ±0.27 
0.41 ±0.11 ±0.03 


13.38 ±2.51 ±0.89 
2.70 ±0.58 ±0.17 
0.44 ±0.10 ±0.03 
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Table 3: The numbers of D* + mesons found with 2 GeV/c < p t D * + < 12GeV/c in bins 
of |r/ D * + | for the three decay modes after background subtraction. The efficiency is listed 
separately for direct and single-resolved processes. The differential cross section in bins of 
the pseudorapidity |^ D * + | of the D* + for each considered D* + decay mode is given together 
with statistical and systematic errors. 



\r] range 


iV found 




D*+ -► (K-vr+)7r+ 


D*+ -► (K-vr+vr )7r + 


D*+ -► (K-vr+7r-7r+)vr+ 


0.0-0.5 


49.2 ±8.9 


21.8±6.8 


51.1 ±10.0 


0.5-1.0 


50.8 ±8.3 


26.4 ±7.6 


45.8 ±9.5 


1.0-1.5 


56.4 ±7.9 


18.5 ±6.3 


29.3 ±7.6 





Efficiency for direct process ef"„ + (%) 


D*+ -► (K-vr + )7r+ 


D*+ -► (K-vr+7r u )7r + 


D*+ -► (K-vr + 7r-7r + )vr + 


0.0-0.5 
0.5-1.0 
1.0-1.5 


41.71 ±0.19 
39.07±0.19 
27.88 ±0.17 


5.45 ±0.09 
5.24 ±0.08 
3.72 ±0.07 


20.90 ±0.16 
19.70±0.16 
12.19 ±0.13 








Efficiency for single-resolved process e|^n*+i(%) 


D*+ -► (K-7r+)vr + 


D*+ -► (K-vr+7r°)7r + 


D*+ -► (K-vr+7r-7r+)vr+ 


0.0-0.5 
0.5-1.0 
1.0-1.5 


37.55 ±0.20 
34.74 ±0.19 
24.08 ±0.16 


4.53 ±0.08 
4.16±0.08 
2.75 ±0.06 


17.31 ±0.16 
15.95±0.15 
9.66±0.11 








da/d\r] B ^\ [pb] 


D*+ -► (K-vr+)7r+ 


D*+ -► (K-vr+vr )7r + 


D*+ -► (K-vr+7r-7r+)vr+ 


0.0-0.5 
0.5-1.0 
1.0-1.5 


13.33 ±2.40 ±0.85 
14.78 ±2.40 ±0.86 
23.22 ±3.24 ±2.10 


12.86 ±4.02 ±1.10 
16.48 ±4.75 ±1.36 
16.59 ±5.63 ±1.73 


14.40 ±2.80 ±1.00 
13.81 ±2.87 ±0.91 
14.35 ±3.70 ±1.31 



Table 4: The combined differential cross sections, da/dp} 3 * and d<j/d|?7 





range [GeV/c] 


2-3 


3-5 


5-12 


dcr/dp^ [pb/(GeV/c)] 


13.50 ±1.51 ±1.01 


3.61 ±0.34 ±0.21 


0.32 ±0.04 ±0.02 






D*+ 1 

rj range 


0.0-0.5 


0.5-1.0 


1.0-1.5 


dcr/d|77 D * + | [pb] 


13.62 ±1.65 ±0.94 


14.65 ±1.71 ±0.94 


18.93 ±2.23 ±1.75 
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Table 5: Sources of systematic uncertainty on the differential cross sections. 



Source 


Estimated uncertainty 




1 U.U U . L ± J /o 


K/7T selection 


(0.5-5.7)% 


Accepted mass range for D° 


< 0.16%, neglected 


D* + selection 


(0.8-2.1)% 


D* + from annihilation events 


< 1%, neglected 


bb background subtraction 


(1.2-3.4)% 


Fraction of direct /resolved r,ar/r res 


(0.3-3.4)% 


BR(D*+ -> D°7T + ) 


2.0% 


BR(D° -> K-7T+) 


2.3% 


BR(D° -> K-tt+tt ) 


6.5% 


BR(D° -> K-7T+7r-7r+) 


5.3% 


Statistical limitation in Monte Carlo 


(0.5-2.3)% 



Table 6: The numbers of D* + mesons found in the acceptance range 2GeV/c < pf* < 
12GeV/c and |^ D * + | < 1.5 for the three decay modes after background subtraction. The 
efficiency is listed separately for direct and single-resolved process. The visible cross 
section a^* + for each considered D* + decay mode is given together with statistical and 
systematic errors. 





(K-7T+)7T+ 


(K-7T+7r )7T + 


(K 7r + 7r 


iV found 


156.4 ± 14.9 


67.4 ± 12.3 


128.4 ± 16.3 


6dir(%) 

e res (%) 


36.47± 0.1 
31.68 ± 0.1 


4.81 ± 0.05 
3.76 ± 0.04 


17.71 ± 0.09 
14.07 ± 0.08 


(Pb) 


24.68 ±2.35 ±1.47 


23.04 ±4.21 ±1.91 


21.76 ±2.76 ±1.41 
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a) Direct process 



b) Single-resolved process 



Figure 1: Main contributions to charm production in 77 events. 



ALEPH 




Am = m D *+ - m D [GeV/c ] 



Figure 2: Mass difference of reconstructed D* + and D° candidates for all considered 
D° decay modes together. The points show data, the error bars represent statistical 
uncertainties, and the solid curve indicates the result of an unbinned maximum likelihood 
fit. 
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Figure 3: Mass difference of reconstructed D* + and D° candidates for three considered 
D° decay modes separately. The points show data, the error bars represent statistical 
uncertainties, and the solid curves indicate the result of unbinned likelihood fits. 
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Figure 4: /Wyis distribution for reconstructed D* + events. The points with error bars 
show data. The relative contributions from direct (shaded histogram) and single-resolved 
(open histogram) processes are extracted by means of a fit. 
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ALEPH 




g§g RES NLO (Massless approach) 
PDF: GRV and AFG 
m c = 1.5 GeV, (J. R =(J. F /2=m T /2 and 2m T 
Default setup: GRV, m c = 1.5 GeV, (i R =(ip/2=m T 



- X//A FO NLO (Massive approach) \ 
PDF: GRS 

m =1.2 GeV and 1.8 GeV 

c 

|J. R =2m T for direct, and nij/2 for resolved 
Default setup: GRS, m c = 1.5 GeV, (J. R =(J.p/2=m T 

i I i I i I i I i I i I i I i I i I i 

2 3 4 5 6 7 8 9 10 11 12 

|^ + [GeV/c] 

Figure 5: Differential cross section da/dp°* + for the inclusive D* + production. The points 
show the combined differential cross sections from the three decay modes under studies. 
The error bars correspond to the quadratic sum of statistical and systematic uncertainties. 
The data are compared to the fixed-order (FO) NLO [12 and the resummed (RES) NLO 
[T3] calculations shown as the solid and dashed lines, respectively. The shaded bands 
represent the theoretical uncertainties of these calculations. 
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Figure 6: Differential cross section da/d|?7 D | for the inclusive D* + production. The 
points show the combined differential cross sections from the three decay modes under 
studies. The error bars correspond to the quadratic sum of statistical and systematic 
uncertainties. The data are compared to the fixed-order (FO) NLO and the resummed 
(RES) NLO [13] calculations shown as the solid and dashed lines, respectively. The shaded 
bands represent the theoretical uncertainties of these calculations. 
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Figure 7: The total cross section for charm production a(e + e~ — > e + e~cc) versus the 
centre-of-mass energy yfs of the e + e~ system. The measurement of this analysis is shown 
as a square. The band represents the NLO QCD calculation [1. The results obtained by 
L3 and OPAL using D* + are represented in [21] and [201 > respectively. The L3 measure- 
ments using lepton tag can be found in [221- The values for TASSO, TPC/27, JADE, 
AMY, and VENUS are taken from [23]. 
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